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Abstract 

Gravitinos and hidden sector fields often cause a cosmological disaster in supersymmetric 
models. We find that a model with gravitational gauge mediation solves such a problem quite 
naturally. The /i-problem is also absent in the model. Moreover, the abundance of gravitinos 
explains correct amount of dark matter of the universe. The dark matter abundance can be 
calculated without detailed information on the thermal history of the universe such as the 
reheating temperature after inflation. 



1 Introduction 



Supersymmetric standard model is advertised as a great successful theory in particle physics 
and cosmology. The hierarchy problem is completely solved by the presence of superpartners. A 
perturbative nature of the model allows us to calculate gauge coupling evolutions in high energy 
and remarkably the three coupling constants meet at very high energy. This is a strong indication 
of supersymmetric grand unification In cosmology, the neutralino, the superpartner of the 
gauge or Higgs fields, can account for dark matter of the universe. The thermal relic abundance 
of the neutralino is in the right ballpark for masses of order 100 GeV. 

However, when we look at actual models of supersymmetry breaking and its mediation, we 
encounter many unnatural aspects or even inconsistencies. In gravity mediated supersymmetry 
breaking scenario there is a serious cosmological difficulty |^ ^ as well as problems of 
flavor changing neutral current (FCNC) and CP violating processes. A supersymmetry breaking 
sector necessarily requires to contain a gauge singlet field whose vacuum expectation value in 
the F-component gives masses to gauginos. After inflation, the scalar component of the singlet 
field starts coherent oscillation. Once it dominates over the energy density of the universe, it 
either spoils the success of Big-Bang Nucleosynthesis (BBN) or overproduces gravitinos H] 
(see 1^, for earlier works). Therefore, for the mechanism of thermal relic neutralinos as dark 
matter to work, we need extra assumptions on the cosmological history so that we can avoid the 
oscillation to dominate the energy density. Models with anomaly mediation ^ can prevent 
such oscillation since singlet fields are not necessary to generate the gaugino masses. However, 
the annihilation cross section of the lightest neutralino, the Wino in this scenario, is too large 
for the thermal relic to naturally explain the dark matter abundance. There seems to be no 
coherent picture for cosmology of supersymmetric models. 



Models of gauge mediation |10| , |11] , |12(| also suffer from serious difficulties. The neutralino 
cannot be a stable particle since the gravitino is always the lightest supersymmetric particle 
(LSP) in this scenario. Moreover, enhanced production cross sections of gravitinos bring up 



a danger of overproduction in the thermal bath 14, |1^, |16|]. This restricts the reheating 
temperature of the universe to be low enough, which makes baryogenesis difficult. There 
is also a naturalness issue in the electroweak sector, i.e., the //-problem. We need to give 
a supersymmetric mass term, //-term, for the Higgs fields which needs to be similar size to 
supersymmetry breaking terms. Unlike in gravity mediation models, we cannot relate this 
/i-term to the gravitino mass via the Giudice-Masiero mechanism |[T7[| because the gravitino 
mass is too small. It is possible to assume a direct coupling between the Higgs sector and a 
supersymmetry breaking sector to obtain the /t-term, but it often induces too large B-term. [O. 
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We cannot say that the minimal supersymmetric standard model (MSSM) is a successful 
model without having a consistent picture for supersymmetry breaking/mediation. It has been 
more than twenty years since discussions of supersymmetry phenomenology started, and until 
recently it has served as the most promising candidate of physics beyond the standard model. 
However, recent various progress in astrophysical observations and theoretical calculations such 



as the measurement of the matter density of the universe |18] and BBN constraints on the 
gravitino abundance put a threat to many supersymmetric models. At the same time, these 
are new strong hints which guide us to the true theory. In this paper, we find that a simple 
model of gauge mediation actually provides us with a viable and natural scenario for cosmology 
and phenomenology. 

The model we discuss consists of a singlet field S and messenger fields q and q. The singlet 
field S breaks supersymmetry spontaneously and that is mediated to the MSSM sector by 
loop diagrams of the messenger fields. The supersymmetry breaking vacuum is stabilized by 



gravitational effects |2C]. This is the entire model. Under reasonable assumptions, the S 
field starts coherent oscillation in the early universe. If we ignore interactions of S with the 
messenger fields, decays of the S field either destroy BBN or overproduce gravitinos, which is 
the same problem in gravity mediation scenarios. A renormalizable coupling to the messengers, 
however, drastically modifies the decay property of S. Especially, S can decay into the MSSM 
particles well before BBN, and, moreover, the branching fraction to the gravitinos is significantly 
suppressed. With gaugino masses of O(IOO) GeV, gravitinos produced by the S decay can be 
dark matter of the universe when the gravitino mass is = O{10) MeV or 0(500) MeV, which 
corresponds to the mass of S to be nis = O(IOO) GeV and 0(400) GeV, respectively. In fact, 
this is the most interesting parameter range from the consideration of the /i-problem. The fi- 
problem can be solved by the mechanism of /i-term driven supersymmetry breaking scenario [^] 
when ms ~ 0{fi). 

In the next section, we explain the gauge- mediation model we consider. The concept of 
/i-term driven supersymmetry breaking is also discussed. In Section ^, we discuss cosmological 
evolutions of the S field. We show that S rolls into a preferable (supersymmetry breaking) 
vacuum in generic situations. The decay property of S is discussed in Section ^. The BBN con- 
straint on parameters and the abundance of non-thermally generated gravitino dark matter are 
calculated. The discussion on the initial condition for the S oscillation is given in Appendix 
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2 A simple model of gauge mediation 



2.1 Gravitational gauge mediation 

The model we study has a very simple structure. We simply introduce messenger fields to any 
supersymmetry breaking model of the O'Raifeartaigh type, i.e., models in which the super- 
potential contains a linear term of a gauge-singlet chiral superfield. The low-energy effective 
Lagrangian is defined by the following Kahler potential and superpotential: 

K = qU + qk+S^S-^-^^ + --- (1) 

W = m^S - XSqq + c , (2) 

where 5 is a gauge singlet field in the O'Raifeartaigh model and q and q are the messenger fields 
which have quantum numbers of the standard model gauge group. We introduced an interaction 
term XSqq where A is a coupling constant. The constant term c contributes to the cosmological 
constant through — 3|VKp term in a scalar potential. The non-renormalizable term — (S'^S')^/A^ 
in the Kahler potential represents the effect of the massive fields (in the O'Raifeartaigh model) 
which have been integrated out. Although this model does not break supersymmetry in the 
limit of the Planck scale Mpi to infinity, a supersymmetry breaking vacuum appears far away 
from the origin of S after taking into account supergravity effects [^]. At the supersymmetry 
breaking vacuum, 

and 

Fs = m^ . (4) 

The messenger particles obtain masses, niq = X{S), in this vacuum. This vacuum is meta-stable 
when 

m^Ml 47r A , , 

where Nq is the number of components in q, e.g., Nq = 5 for a pair of q and g in 5 and 5 
representation of SU(5). 

The constant term c in the superpotential is chosen so that the cosmological constant is 
cancelled at the supersymmetry breaking vacuum, i.e.. 



c 



m3/2M|i , (6) 
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where the gravitino mass is given by 



2 

m 



The scalar component of S obtains a mass by supersymmetry breaking effects: 

= ■ (8) 

Supersymmetry breaking by Fs = induces gaugino masses by one-loop diagrams with 
messengers q and q. It is calculated to be: 

_ g^N Fs _ g^N 2^/3m^Mpl 
""^ (47r)2(5) (47r)2 A2 ' 

where is the number of the messenger particles (A^ = 1 for a pair of 5 and 5 representation 
of SU(5)). The phenomenological requirement that the gaugino mass should be O(IOO) GeV 
determines the ratio m/A. 

Among three parameters m^, A and A in this model, low energy physics is not very sensitive 
to A. It only appears through logm-g where rrig = X{S) is the messenger scale of gauge mediation. 
With fixed gaugino masses, we are left with a single parameter A. 

2.2 /Lt-term driven supersymmetry breaking 

There is a particularly interesting scale for A from the consideration of UV completions of the 
theory. Among various possible UV completions (see Ref. ||2^ for various possibilities), the 
simplest one is to identify the field S with a meson field of a strongly coupled gauge theory. For 
example, in an SO (5) gauge theory with a quark T and a superpotential term: 

W = uTT , (10) 
the low energy effective theory below the dynamical scale A is simply 

W = nAS , (11) 

in one of two branches of the theory. Supersymmetry is broken by Fs = /xA. This vacuum is 



found to be meta-stable recently by Intriligator, Seiberg and Shih |22]. The Kahler potential is 
unknown but can be parametrized to be 

K = StS-2(^ + ...^ (12) 
A 

If T] is positive, this is exactly the same effective theory we discussed in the previous subsection. 
The scale A is identified with kjr^l'^ . The messenger fields couple to S through W 9 TTqqjMt^ 
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in the superpotential with the coefficient of the operator. The A parameter in Eq. (^) is 

then A = A/M*. 

In this model, the /i-term in Eq. ( [lO| ) can actuahy be the ^u-term for the Higgs fields W B 
fiHuHfi in models with composite Higgs fields |21|. The Higgs fields and are meson fields 



Hu ~ (QT) and Hd ~ {QT) where Q and Q carry the standard model quantum numbers. The 
mass term for T in Eq. ( [lO| ) therefore induces a mass term for mesons. (For detailed discussion, 
see ||2l|.) Since S and the Higgs fields Hu and Hd are strongly coupled, we expect to have terms 
in the Kahler potential: 

K 3 ^S^SHIHu + ^S^SH\Hd . (13) 
A A 

With Fg = yuA, these terms induce the soft supersymmetry breaking terms 

~ "^L ~ /^^ > (14) 

independent of the dynamical scale A. This provides a solution to the /x-problem in gauge 
mediation. The /i-term is related to the supersymmetry breaking parameters because it is the 
source of the supersymmetry breaking. 

A combination of the gravitational gauge mediation and /i-term driven supersymmetry 
breaking provides a completely natural scenario for gauge mediation. The m? parameter in 
Eq. (|2|) is 0{^A) in this scenario (under a reasonable assumption that rj = 0(1)), and it fixes 
the dynamical scale A(~ A) by the requirement /i ~ nix ~ O(IOO) GeV from Eq. For 
having this relation, the factor Mpy/A needs to compensate the one-loop factor g'^ N / (Att)'^ , 
which determines the dynamical scale to be of the order of the GUT scale A ^ 10^^ GeV . This 
is indicating a close connection between dynamics of supersymmetry breaking and spontaneous 
symmetry breaking in grand unified theories. In fact, the explicit model in Ref. [21| is based 



on a GUT model. A meson field made of Q and Q introduced before plays a role of the Higgs 
field which breaks GUT symmetry dynamically. This mechanism simultaneously provides us 
with a solution to the doublet-triplet splitting problem with sufficiently suppressed dimension 
five operators for proton decays. 

With ~ fiA, the mass of scalar component of S, ms, is 0(/i) ~ 0(100 — 1000) GeV (see 
Eq. (|8|)), and the gravitino mass is 0(0.01 — 1) GeV. We will use ms and ^ ^ee parameters 
instead of m? and A. Relations among them are: 

= (1 X 10^ GeV)2 x ( ""^^^ ] , (15) 
^ ^ V500 MeV/ ^ ^ 

(1.) 
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The Bino mass is expressed in terms of these two by 

"l3/2 



m . = 200 GeV x N ( ) ^ (^7) 

^ V 400 GeV/ VSOOMeV/ ^ ^ 

There are no 0(1) ambiguities in these relations. 



3 Cosmological evolution of S 

As we have seen in the previous section, the S field is stabilized at the supersymmetry breaking 
vacuum by supergravity effects. However, there is a true and unwanted supersymmetric vacuum 
at S* = and q = q = y^m^/A. It is thus non-trivial whether S can settle to the supersymmetry 
breaking minimum in the thermal history of the universe. This question is addressed in Ref. [^] 
where it is concluded that thermal effects make the supersymmetric vacuum more attractive 
once the messenger particles are thermalized. In our scenario, it is possible that the temperature 
of the universe never exceeds the mass of messenger particles so that thermal effects are not 
important. 

We show in this section that S runs into the supersymmetry breaking minimum under 
reasonable assumptions. It is found to be rather difficult to fall into the supersymmetric 
minimum S = because it can be thought of a small hole in the vast S space. 

The relevant scalar potential of the S field in the IS"! < A region is given by 

V{S) = m' + ^ log +•••)- {2V^ml,Mp,S + h.c.) . (18) 

The potential for IS"! > A region depends on a UV completion, but it is not important for 
our discussion. The logarithmic term is induced at one-loop level due to the interaction with 
messenger particles. The supersymmetry breaking vacuum given in Eq. (^) is determined by a 
balance between the mass term and the linear term. The logarithmic term is important near 
S = 0. The S field, however, feels much stronger force toward the supersymmetry breaking 
vacuum at most points on the complex S plane since it is supported by a quadratic potential. 
This is why the supersymmetry breaking vacuum is more attractive at a low temperature. If 
the value of S is not at the origin after inflation, S rolls down to the S ^ vacuum. 

The minimum of the scalar potential during inflation is not at S" = or S = {S), but is likely 
to be around l^l ~ A (see Appendix^ for detailed discussion). Once inflation starts, S is trapped 
at the minimum 5*0 = 0(A) and freezes until the end of inflation. After inflation, the value of 
S remains fixed at S ~ 5*0 for a while until the mass term ms becomes important compared to 
the Hubble mass term of 0{H), and eventually, it starts to oscillate about the supersymmetry 
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Figure 1: A typical evolution of the S field (left) and time dependence of the amplitude |5| 
(right) for ms = 400 GeV, mg/s = 500 MeV and A = 10^^. 

breaking minimum. In general, the initial value Sq is not aligned to the {S) direction, which 
makes motion of S almost confined in a straight line on the complex plane. Therefore, S never 
approaches to the supersymmetric minimum. The amplitude of the oscillation gradually decays 
due to the cosmic expansion and S settles to the supersymmetry breaking minimum. 

In Fig. |l|, we show a typical evolution of the S field for ms = 400 GeV, 171^/2 = 500 MeV 
and A = 10~^. Here, the initial condition we took is Re(S') = Im(S') = A. As we expected S 
oscillates along an almost straight line about the supersymmetry breaking minimum (see also a 
closeup in the left panel). We also show time dependence of the amplitude |S| during oscillation 
in the right panel. The figure shows that the impact parameter to S" = during the oscillation 
is of 0{{S)) and the amplitude of the oscillation decays gradually. Eventually, S settles to the 
supersymmetry breaking minimum. 

In the above analysis, effects of the logarithmic potential in Eq. ( p!8| ) are not important as 
long as A < 10"^^"^^. It does attract 5 to the supersymmetric minimum, but the attractive 
force is only important when S approaches very close to the origin. With an impact parameter 
of 0{{S)), the effect is negligible. 

The messenger fields have masses, niq = XS, when S is away from the origin. With the large 
impact parameter of the oscillation of S, q and q are always stabilized at q = q = 0. 

Although they are always massive, time variation of the masses by the S oscillation causes a 
non-perturbative production of the messenger fields |2^. Once they are copiously produced, it 
gives rise to a force to attract S to origin so that the energy density is reduced by minimizing the 
mass of q However, since the variation of niq is quite adiabatic (i.e., rhq ~ msiriq <C rUq), 

the effect is exponentially suppressed. 

*We thank J. Wacker for discussion on this point. 
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We need to check here that the messenger particles do not enter thermal equilibrium for 
the above discussion to hold. If they are produced in the thermal bath, it induces a thermal 



mass term to S which attracts to 5 = [23|. However, as long as temperature is always 
(including the inflaton oscillation era) lower than niq = XS, production of messenger particles is 
also exponentially suppressed. With the large impact parameter {0{{S))) and initial amplitude 
(0(A)) of this scenario, the condition can be satisfied in most inflation models. 

4 Gravitino dark matter 

We now follow the cosmological evolution after the S oscillation. A brief summary of the 
history is the following. Because of a long lifetime of the S field, its oscillation energy density 
eventually dominates over the universe unless the reheating temperature after inflation is very 
low. Decays of S reheat the universe to a temperature of 0(1 — 100 MeV). The large entropy 
production from the S decay dilutes preexisting radiation and matter, including thermally 
produced superparticles, gravitinos and also even baryon asymmetry of the universe. The main 
decay channel of the S field is into two Higgs bosons (5 — > hh) or two gluons {S gg) depending 
on ms- Radiation produced by the S decay takes over the energy density of the universe below 
this decay temperature and therefore standard BBN can happen after temperature drops down 
to O(MeV). Decays of S produce gravitinos through a rare decay mode S — > The 
produced amount is remarkably consistent with the observed dark matter density of the universe 
in the parameter regions of our interest. 

We start discussion with deriving the interaction Lagrangian of S, with which we can 
calculate the decay width of S and the branching ratio of the gravitino mode. We will then 
discuss the gravitino abundance of the universe. 

4.1 Interaction Lagrangian of the S field 

The S field couples to the MSSM particles only through loop diagrams with messenger particles 
q and q. The effect is encoded in the S dependence of low energy parameters. 

The low energy values of the gauge coupling constants depend on S through log S since it 
changes the scale at which q and q are integrated out. The kinetic term of the gauge bosons are 

Ain = -^^Pf^-P^" + h-C. , (19) 



where 



1 2{bL + N) XS 2bL , fiR 

■ log 7 7:r^ log — , (20) 



g^iS) gi (47r)2 " (47r)2 ^ XS 
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with go being the gauge couphng constant at a scale Ah- The factor 6^ is the beta- function 
coefficient below the messenger scale and /i_R is the renormalization scale. After canonically 
normalizing the kinetic terms of gauge bosons, we obtain the interaction terms from Eq. (|l9|): 

g'^N 1 1 



(47r)2(5) 4 



SF^,F>"' + h.c. 



21 



The interaction terms with gauginos can be obtained in a similar way. From the S dependence 
of the gaugino mass terms: 



augino 



^mx{S)XX + h.c. , 



with 



we obtain 



mxiS) 



g'^N m? 



(22) 



(23) 



(24) 



Comparing with Eq. (|2lD, this gives larger contribution to the decay width of 5 by a one-loop 
factor when ms ~ ttt-a- 



The same conclusion can be obtained for scalar fields. The soft mass terms are 



>Cscalar = -mj{S)Pf , 



with 



m'jiS) 



(47r)^ 



2C2N 



m 



S 



C2 is the quadratic Casimir factor. We can read off interaction terms as follows: 



2 



-f 



(25) 



(26) 



(27) 



Again, this gives larger contribution by a one-loop factor compared to Eq. (^T]). For the Higgs 
fields, the scalar masses have two origins, i.e., one from the gauge mediation and the other 
from Eq. (^). For the calculation of couplings between S and the Higgs fields, we should use 
the gauge- mediation contribution, i.e., Eq. (^), for in Eq. (^). It is amusing that the S 
field dominantly decays into supersymmetric particles and the Higgs fields if it is kinematically 
allowed. 

Coupling to gravitinos can be obtained by simply looking at the Lagrangian. The largest 
contribution for the S decay comes from the coupling to the longitudinal mode of the gravitino 
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(the fermionic component of S, s). From the — (S'^S)^/A^ term in the Kahler potential, we 
obtain to ss coupling such as 

A/2 = -^S^s's + h.c. ^ -i!^5tv;3/2A/2 + h.c. (28) 

Therefore, the decay width is suppressed by 0{{m^/2/mx)'^) compared to the gaugino/scalar 
modes. Also, there is a suppression of 0(^,3/2/^,5) compared to the gauge boson mode. This 
effect is important for the gravitino abundance. 

4.2 S decays and gravitino production 

In the following we consider two cases (A): ms > 2mh, where S can decay into two Higgs bosons 
and (B): ms < 2mfi, where S dominantly decays into two gluons. We find that in both cases 
gravitinos from the S decay in early time naturally explain the dark matter component of the 
universe. 



Case (A): The S — > /i/i decay is open {ms > ^m^). 

We further assume here that S dominantly decays into two Higgs bosons. This assumption is 
valid as long as ms/2 is smaller than other 811(2)/, or SU(3)c charged superparticle masses 
(except for the Higgsino). This will be justified later. From the experimental lower bound on 
the Higgs boson mass m^ > 114 GeV |26|, we obtain ms > 230 GeV. The decay width of S in 
this case is given by 

_ xjjN'^ m| / ms V 
^ ~ 15367r V"i3/2/ ' 

where xh is a two-loop factor: 

g2 and gy are the gauge coupling constants of SU(2)i and U(l)y gauge interactions. The lifetime 
of S is then 

rs = 5x 10-5 sec x N'^ f-^V'' (^^Y • (31) 



400 GeV/ V500 MeV 

With this relatively long lifetime and the large initial amplitude of S, it is reasonable that 
the coherent oscillation of S dominates over the energy density of the universe before S starts 
to decay. Although it is possible that S domination does not happen by assuming a presence of 
another long-lived matter density (such as the inflaton oscillation), we do not consider such a 
case here. 
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The decay of S produces radiation and it reheats the universe. The reheating temperature 
after the S decay is calculated to be 

. 0.5 X vfw;; = 90 MeV x A. (j^)'"' (siiiw)" ' (^^> 



This temperature needs to be larger than 2 MeV in order for standard BBN to happen |27|. The 
baryon asymmetry and the dark matter component of the universe also needs to be generated 
before BBN. Since primordial radiation and matter are significantly diluted by the entropy 
production from the S decays, it is non-trivial whether we can consistently obtain those two 
components. We explain here the dark matter (gravitino) production from the S decay and we 
will discuss a possible mechanism for baryogenesis later. 

Under the assumption that the S'-domination happens, the number density of the gravitinos, 
723/2, from S decays is given by 

^ = |^i^3/. X 2 , (33) 
s 4 ms 

where is the branching fraction into two gravitinos and s is the entropy density of the 
universe. Note here that the number density can be calculated only with the two parameters in 



(34) 



the Lagrangian, ms and The partial decay width is calculated from Eq. (|2 



^Pl V"i3/2, 

Therefore, the branching fraction i?3/2(= '^3/2/'^h) is given by 

By. = 2 X 10-6 X N-^ ( y ( "^"^ ]" . (35) 

V400 GeV/ V 500 MeV/ ^ ^ 

The energy density is then 

Comparing with the observed dark matter density of the universe, ncDuh^ = 0.10 ± 0.02 
(2(t) , we find that the gravitino is a perfect candidate of dark matter in this scenario. 

If S can decay into two Binos and if the Bino is the next to lightest supersymmetric particle 
(NLSP), the produced Binos later decay into gravitinos. This possibility is excluded by the BBN 



constraints on the Bino decay [19|. The conclusion is the same for the case where the Higgsino 
is lighter than the Bino. This justifies the assumption that S mainly decays into two Higgs 
bosons. Other SU(2)/, or SU(3)c charged supersymmetric particles are heavier than the Bino 
in gauge mediation. 
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In fact, if the scalar tau lepton (stau) is lighter than the Bino, which is possible if there is a 
large ^f-term from the direct coupling between the Higgs fields and S in the Kahler potential, 
a small portion of the ms > 2m ^ region becomes viable. Staus can reduce the number density 
by the pair annihilation process after the production from S decays. The energy density of the 
stau before its decay is approximately 

n^h^^nfh^'^, (37) 

where Tj is the freezing out temperature, Tj ~ mf/20, and f^NLSP^^(~ 0.002(mf/100 GeV)^) 
is the energy density of the stau obtained by the standard calculation of the thermal relic 
abundance. In the region of our interest, where the gravitinos from the S decay explain dark 
matter of the universe, the above amount is on the border of the BBN constraints. Since the 
viable region appears only in a special circumstance and also the region is quite small, we do 
not consider this possibility in the rest of the paper. 



Case (B): The S ^ hh decay is closed {ms < ^m^). 

In this case, the dominant decay channel is into two gluons, S gg. The decay rate is calculated 
to be 



^ 967r \m^/2 



with a one-loop factor 



- (47r)2 

From this, the reheating temperature after the S decay is given by 



(38) 



fl2 

Xo = = 9.4 X 10"^ . (39) 



VTOGeV/ VlSMeV/ ^ ^ 



The branching fraction of the gravitino mode is 

B,,, = 5xW-'xN-'(-^] {—^] . (41) 



By using the same formula in Eq. (^), we obtain 

n,,,h' = 0.1 X N-^ f^I^y^' (J!^) . (42) 
^'^ V70GeV/ VlSMeV/ ^ ^ 

The dark matter density can be explained by the non-thermally produced gravitinos also in this 
case. 
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Figure 2: The parameter regions of the model {N = 1) which account for the dark matter 
component of the universe. The correct amount of gravitino dark matter can be obtained 
within the allowed range of parameter space. Below the dot-dashed line in the left figure, the 
decay mode of S into two Higgs boson is closed (Case (B)). The constraint from the lower bound 
on the decay temperature Td > 2 MeV is indicated. The contour of = 5 MeV is also shown. 
In the right figure, the decay mode 5 — > /i/i is open above the dot-dashed line (Case (A)). The 
decay mode into two Binos is open in the right-lower region. 

We show in Fig. ^ the region of the parameter space where the non-thermally produced 
gravitino explains the dark matter component of the universe. We took = 1 and used 
m3/2 and as free parameters of this two-parameter model, and showed contours of ms (see 
Eq. (|l^) as dashed lines. The values of 1^3/2/1^ expressed in terms of those two are 



(43) 



^3/2 / \l/4 

.15 MeV/ V200 GeV. 
respectively. The right figure corresponds to the Case (A) where the Higgs boson mode is 
the dominant decay channel (above the dot-dashed line). The constraint from > 2 MeV 
is satisfied in the figure once the Higgs mode is open. The lower bound on the Bino mass 
(m^ > 160 GeV) is obtained from the requirement that the Bino mode should be kinematically 
forbidden. The gravitinos can explain the dark matter density when ms ~ 300 — 500 GeV, 
which is exactly the preferred values from the /x-term driven supersymmetry breaking scenario 
and grand unified theories. 

Another consistent region is found in the left figure, where the gluon mode is the dominant 
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decay channel (Case (B)). The constraint from BBN, > 2 MeV, excludes the lower mass region 
of in this case. The correct amount of dark matter can be obtained with ms ~ 100 GeV. 
This is also a natural size for the ^-term driven supersymmetry breaking scenario. 

For a larger value of N, we obtain a larger lower bound on m^, and the allowed region of 
is iiot very sensitive to for Case (A) and moves to larger values linearly for Case (B). 

We comment on other possible sources of gravitinos. Thermally produced gravitinos |l^, |l^, 
|l5|, are significantly diluted by the entropy production from S decays. It is diluted by a factor 
of Tfi/T^om where Tdom the temperature at which S oscillation starts to dominate (see discussion 



in the Subsection 4.4). If S oscillation starts during the reheating process after inflation, which 
is the case in our scenario in order not to thermalize the messenger fields, the amount of thermal 
gravitinos does not depend on the reheating temperature after inflation. It is given in terms of 
the parameters of the model and the initial amplitude |5o| by 
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This is sub-dominant when |5o| ~ A. 
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(44) 



Thermal relic of NLSPs also generates gravitinos through their decays |28]. Again entropy 
production reduces the amount of NLSPs by a factor of Max[(T(i/Tj)^, T^^/ (Tkom^/)^^^) ?d/Tdom]) 
where T/(~ mNLSp/20) is the freezing out temperature of pair annihilation processes [^9|. The 
super WIMP component of gravitino dark matter is therefore negligible in the parameter regions 
of our interest. 



4.3 Are gravitinos cold? 

Since gravitinos are non-thermally produced through the decay of a heavy particle S, they are 
relativistic when produced, but they subsequently slow down as the universe expands. The 
gravitinos need to be sufficiently non-relativistic before the structure formation to start. The 
free-streaming scale of the dark matter is constrained to be less than 0(1) Mpc by the matter 
power spectrum inferred from Lyman-a forest 



We estimate the free-streaming scale Aps of the gravitinos. We find that it is sufficiently 
small, but interestingly it may be large enough to be distinguished from that of cold dark matter. 
The free-streaming scale at the time of matter-radiation equality icq is defined by 



a{t) 



dt ~ 



ms 



an 
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(45) 
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where Oeq and ao are scale factors at matter-radiation equality and at the time of S decays, 
respectively. Hq is the Hubble parameter at the current epoch. Up to logarithmic dependence, 
we obtain 

"'^'(^7^^^)' (Case (A)), 

(46) 



Afs 



0.004 Mpc X N 



-1 



400 GeV 



For Case (B), Aps is in an interesting range for the small-scale structure problem in ACDM 
model. Although not conclusive, it has been reported that there are discrepancies between 
simulations of structure formation in the cold dark matter scenario and observation at a small 
scale. The observed number of galactic satellites in the inner region (O(Mpc)) of halo is an order 
of magnitude smaller than results of simulations pl| . Also, simulation of the structure formation 



predicts a cuspy profile in the central region of halo [ 32 1 whereas observations such as studies of 
gravitational lensing 



and rotation curves |34] indicate nearly uniform density core structure. 
It has been studied that dark matter with free-streaming scales of 0(0.1 — 1) Mpc may resolve 
these issues because it can smear out the density fluctuation in small scale (See also 

p7| for recent proposals of particle physics models which realize large Aps via non-thermally 
produced dark matter.) 

A free-streaming scale of 0(10 — 100) kpc means that there is a cut-off in matter power 
spectrum below that scale. It might be possible to probe this small-scale structure in future 

1, and 



experiments by using strong gravitational lensing |38] such as in LSST [39| and SNAP 



also by observations of 21cm fluctuations in high red-shift region [41|. 



4.4 Entropy production and baryogenesis 

In the previous section, we have assumed that the oscillation energy density dominates over the 
universe before it decays. We discuss the amount of entropy production and possible baryogenesis 
scenarios here. 

As we have discussed, the S field starts to oscillate when the Hubble parameter becomes 
smaller than ms- There are two cases which we need to consider separately. If the oscillation 
starts during the reheating process after inflation (the inflaton oscillation era), S'-domination 
happens when: 

where T/j is the reheating temperature after inflation. 
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A different result is obtained if S starts to oscillate in a radiation-dominated era (after 
reheating). The temperature at which S starts to oscillate is given by 

Tosc ^ 0.3 X ^/Mpms =^ 8 X 10^ GeV x {^^^) • (48) 



In this case, ?dom is 



The dilution factor A is then 



' dom 



(50) 



The primordial radiation and matter are diluted by this amount. (S-domination does not 
happen if is larger than unity.) With the low decay temperature T^^, dilution effects are 
quite large. For example, if Tr < Tqsc and /So ~ A ~ MquTi the dilution factor is ~ 
10~^(Tr/10^ GeV)~^. This dilutes unwanted relics such as thermally produced neutralinos to 
a negligible level. Since baryon asymmetry is also diluted, primordial baryon asymmetry needs 
to be larger such as ub/s ~ 10^^. It is a non-trivial task to obtain such a large asymmetry 
under a condition that the reheating temperature cannot be too high in order not to thermally 
produce the messenger particles (min[(TRTosc)^^^, Tqsc] ^ ^{S))- 

Examples of possible scenarios for baryogenesis are the right-handed sneutrino inflation sce- 
nario and scenario with right-handed sneutrino dominated universe |^3|. In those scenarios, 
baryon asymmetry is generated via leptogenesis after the decay of right-handed sneutrinos 
as follows: 

!^.,0,2-0.8,.:0-.(^)(^).n4. (51) 

where 5es is an effective CP violating phase and rn-j^s corresponds to the heaviest neutrino mass, 
and Tr denotes the decay temperature of the right handed sneutrino. After dilution by the S 
decay, it reduces to 

^ ^ (0.2 - 1) X 10-1° X ( ) ( 'f°' ' ( ^^"l sin^eff, (52) 

s ^ ' V 100 MeV y V 1015 GeV y VO-05 eVy ^ ' 

independent of Tr for Tq^c > T^. Therefore, the observed baryon asymmetry n^/s ~ (8.7ib0.3) x 

lO^ii ||l^ can be consistently explained with a slightly suppressed initial amplitude |S'o|. Note 

here that if the initial amplitude Sq is suppressed to this level, thermally produced gravitinos 



can give a non- negligible contribution to the matter density (see Eq. (44)), and also the decay 
of thermally produced NLSPs during BBN may put constraints on species of the NLSP and/or 
Tr. 
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5 Summary 



Taking into account gravity effects in supersymmetry breaking/mediation models opens up a 
possibility of very simple gauge mediation. We have considered the simplest model of gauge 
mediation, in which we have only two parameters. One combination of the parameters needs to 
be fixed so that we obtain the gaugino masses to be O(IOO) GeV, and the other has preferable 
range from a simple solution to the /U-problem and grand unification, i.e., 10 MeV < ?7i3/2 ^ 
1 GeV. 

This is the simplest model of gauge mediation, but amazingly free from problems. We have 
shown that S rolls into the supersymmetry breaking vacuum with generic initial conditions. 
The entropy production by the S decay dilutes gravitinos and neutralinos which are thermally 
(over)produced. Supersymmetric particles do not cause the FCNC or CP problems as the 
gravitino mass is small enough. Of the most amazing is the consistency of the dark-matter 
abundance. We can calculate the abundance of the non-thermally produced gravitinos from S 
decays in terms of the two parameters in the model. Within the parameter region of our interest, 
the calculated amount can account for the dark matter of the universe. 
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A Initial condition of S 

As we mentioned in section ^, the minimum of the scalar potential during inflation is likely to 
be around \S\ ~ A. 

The singlet field S may couple to the inflaton sector in the superpotential, 

W = inflaton (l " +•••), (53) 

where ^ is a coefficient of 0(1), and Winflaton denotes the superpotential of a inflaton sector. 
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The scalar potential of S during inflation is then given by 



V{S) ~ ( \S\^ - ^ + ...)+ 



(54) 



where the first term originates from the Kahler potential in Eq. (1) and the second term from 
the superpotential in Eq. ([5^). In the above expression, we have used the Hubble equation, 

tt2 I -^Inflaton I ^ 
11 — , 

3M2j ' 

assuming that the F-component of the inflaton dominates the energy density during inflation. 
This potential has a minimum at 5*0 = 0{A). 

In principle, the ^ parameter can be much larger than 0(1) depending on a UV completion 
of the theory. For example, in a model of Ref. |^l|, the theory above the scale A is a strongly 
coupled conformal field theory, and the ^ parameter is enhanced to be 0(Mpi/A). If this is 



the case, the second term in Eq. (54) pushes the minimum to the Planck scale. However, the 
effective theory in terms of 5 does not make sense for IS*! > A, and it is more appropriate to 
consider the dynamics in terms of the constituent field T. The constituent field T, in contrast 
to S, may have a positive Hubble mass squared around the symmetry enhancement point at 
r = 0. This Hubble mass term attracts T to the origin, and it prevents S ~ TT from going to 
the Planck scale. Therefore, we can expect the actual minimum during inflation does not exceed 
0{A) even if the coefficient ^ is large. 

If ^ is large, higher order terms of S in Eq. ( [5^ ) can be also important. Even in this case, 
the result of our discussion will not change qualitatively. 
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